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Julia Fang^ and Jean-Luc Margot^'^ 

Draft version November 14, 2011 

ABSTRACT 

In the near-Earth asteroid population, binary and triple systems have been discovered with mutual 
orbits that have significant eccentricities as well as large semi-major axes. All known systems with 
eccentric orbits and all widely-separated primary-satellite pairs have rapidly rotating satellites. Here 
we study processes that can elucidate the origin of these spin-orbital properties. Binary formation 
models based on rotational fissioning can reproduce asynchronous satellites on orbits with high eccen- 
tricities and a wide range of separations, but do not match observed properties. We explore whether 
any evolutionary mechanisms can link the spin and orbital parameters expected from post-fission 
dynamics to those observed today. We investigate four processes: tidal torques, radiative perturba- 
tions (BYORP), close planetary encounters, and Kozai oscillations. We find that a combination of 
post-fission dynamics and tidal evolution can explain nearly all the spin-orbit properties in a sample 
of nine well-characterized near-Earth binaries and triples. The other mechanisms may act but are not 
required to explain the observed data. Lastly, we describe evolutionary pathways between observed 
spin-orbital states including synchronous and circular, asynchronous and circular, and asynchronous 
and eccentric configurations. 

Subject headings: minor planets, asteroids: general - minor planets, asteroids: individual (2000 DP107, 
1999 KW4, 2002 CE26. 2004 DC, 2003 YTl, Didymos, 1991 VH, 2001 SN263, 
1994 CC) 



1. INTRODUCTION 

The study of near-Earth asteroids (NEAs) and main 
belt asteroids (MB As) with satellites can yield important 
information about their fundamental phys ical properties 
as well as their formatio n and evolution (jMerline et al.l 
120021 : IPravec et al1l2006[ ). In these multi-component sys- 
tems, analysis of the relative positions between compo- 
nents can quantify their mutual orbits as well as the to- 
tal mass of the system. In the near-Earth population, 
binary and triple asteroids are typically discovered by 
radar (24 out of 37 as of October 2011) during a close 
approach to Earth, which can provide detailed physical 
and orbital informa tion about the system. Studies by 
iMargot et all (|2002D and IPravec et al.l (|2006f ) have deter- 
mined that approximately 15% of NEAs larger than 200 
meters in diameter have satellites. Due to the depen- 
dence in the return signal (where d is the distance to the 
target), radar observations have not identified MBA sys- 
tems so far. Small MBA systems have generally been dis- 
covered through light curve observations whereas larger 
MBA systems have typically been discovered with adap- 
tive optics observations. 

Orbital solutions of NEA systems indicate that some 
of their satellites possess unexplained spin-orbital prop- 
erties including asynchronous rotation, eccentric orbits, 
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and wide separations from their primaries. All known 
satellites with semi-major axes larger than 7 primary 
radii arc asynchronously rotating and all known satellites 
with eccentricities greater than 0.05 are asynchronously 
rotating, suggesting that these properties have a common 
origin. Accordingly, we seek an explanation for these ob- 
served spin-orbital characteristics by examining if any 
evolutionary processes can lead to the observed data. 

Previous attempts to investigate the observed prop- 
erties of NEA systems include tidal evolution as a 
mechanism for eccentri city excitation or de-excitation 
(|Tavlor fc Margo'3l2011l and references therein). Other 
studies have described the orbital evolution of small bi- 
nary asteroids by the binary YORP (BYORP) effect 

HCuk fc BurnsllMOSi : [Cui3[2051 . BYORP is caused by 
the asymmetric rc-radiation of light by an irregularly 
shaped secondary in synchronous rotation with its pri- 
mary, and this effect can cause orbital migration and 
an incr ease or decrease of th e mutual orbit's ccccn- 
tricitv (IGoldreich fc SanI [2001 ICuk fc Nesvorn v 20 ]J; 
IMcMahon fc Scheereal ' IToiOaU d: ISteinberg fc Saril l201]F . 
An alternative process to modify an orbit's semi-major 
axis and eccentricity is through close scattering encoun- 
ters by terre strial planets, as desc ribed in a compan- 
ion paper bylFang fc MargoB (|2011t ) for binaries and by 
iFang et aLl (|2011[ ) for triple svstems. Another possibility 
is that NEA binaries can h ave their ecce ntricities excited 
through Kozai oscillations (jKozailll962[ ). 

In this paper, we examine all of these proposed evo- 
lutionary processes to find a coherent model that can 
explain the observed spin-orbital characteristics of satel- 



120071) . who use the term "asynchronous binaries" for binaries with 
spin-orbit synchronization. If generalization to systems with more 
than one satellite is needed, we affix the terms synchronous and 
asynchronous to the satellites being considered. 



Table 1 

Well-Characterized Near-Earth Binaries and Triples 
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(185851) 2000 DPIOT^ 
(66391) 1999 KW4^ 
(276049) 2002 CE26^ 
*2004 DC^ 
*(164121) 2003 YTr 
(65803) Didymos^ 
*(35107) 1991 VH^ 
(153591) 2001 SN263 #1^;. 
*(153591) 2001 SN263 #2^ 
(136617) 1994 CC #r 
*(136617) 1994 CC #2^ 

Well-characterized NEA binaries and triples as defined in this paper have some known physical and orbital properties: approximate 
component sizes (primary radius secondary radius R.s), masses (primary mass Mp, secondary mass Ms), semi-major axis a, and 
eccentricity e. Known asynchronous satellites are marked with a *. This table shows nominal values adopted for this study as well 
as plausible l—cr uncertainties for a and e. Uncertainties for sizes are r ough ly ~20% and for masses are roughly ~10%. Methods for 
obtaining parameters and uncertainties are described in the text (Section ll.lll . Also shown here is close approach (CA; < 0.1 AU) data, 
including the most recent approach at the time of radar observation and the current MOID (see text) with Earth, given by the JPL 
Small-Body Database. For the triple systems, 2001 SN263 and 1994 CC, we list the inner satellite first and the outer satellite second. 

'IMargot et al.l (f200l) 

fOstro et al. (2006) 

iShepard et al. (2003) 

<^ Taylor et al. (2008) 

•tNolan et al. (2004) 

|Benner et al. (201(|) 

SMargot et al. ( 2008); IPr avec et al.| l|200i) 
''Nolan et al. (2008); Fang et al. (2011) 
IBrozovic et al.. (,2011. 1. .Fang et al.. 12011i 'l 

lites in NEA systems. In the remainder of this section, 
the current population of well-characterized NEA bina- 
ries and triples is presented in Section 11.11 relevant life- 
times are defined in Section II. 2[ and binary and triple 
formation is introduced in Section 11.31 Then, we ex- 
amine each main evolutionary process in turn: Section 
[2] discusses tidal evolution timescales and critical semi- 
major axes pertinent for satellite spin synchronization, 
Section [3] summarizes current theories on BYORP evo- 
lution and their applicability, Section|4]evaluates if plane- 
tary encounters can explain the eccentric and wide orbits 
of asynchronous satellites, and Section [S] discusses Kozai 
resonance and its timescales. We find that tidal evolu- 
tion can explain the observed spin-orbital characteristics 
of nearly all NEA systems, and we discuss possible evolu- 
tionary pathways between observed spin-orbital states in 
Section [6l Wc summarize this study and its implications 
in Section [71 

1.1. Sample of Well- Characterized Binaries and Triples 

We compile a sample of well-characterized NEA sys- 
tems (Table [1]) that consists of 7 binaries and 2 triples. 
We will refer to this sample throughout this paper as we 
study the spin-orbital origin of these systems. This sam- 
ple includes all NEA systems with known estimates of 
system mass, semi-major axis, eccentricity, and compo- 
nent sizes. In practice, only systems observed with radar 
fall in this class. 

Primary and secondary sizes of these NEAs are ob- 
tained from radar shape modeling (when available), or 
from range extents estimated from radar images. System 
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masses are derived from orbital solutions, which are com- 
puted based on measurements of range and Doppler sep- 
arations. For triple systems, individual mass estimates 
are obtained through dynamical 3-body orbital fits. For 
some binaries, the masses of the primary and of the sec- 
ondary have been directly estimated from the observed 
reflex motion. For all others, the system mass is appor- 
tioned to the individual components by using size ratios 
and a common density assumption. 

The semi-major axes and eccentricities of these well- 
characterized NEAs are obtained through orbital fits to 
radar astrometry, and we list them in Table [T] and plot 
them in Figure [T] iBrozovic et al.l ()201lD have reported a 
list of asynchronous satellites that are rapidly rotating, 
which includes the following satellites in our sample of 
NEA systems: 2003 YTl, 1991 VH, 2004 DC, and the 
outer satellites of 2001 SN263 and 1994 CC. These asyn- 
chronous satellites are marked with asterisks in the first 
column of Table [1] and represented by unfilled circles in 
Figure [TJ In this sample, asynchronous rotators include 
all satellites with adopted eccentricities greater than 0.05 
and all satellites with semi-major axes greater than 7 pri- 
mary radii. These correlated spin-orbital properties are 
likely due to the decreasing effects of tidal dissipation 
(which can synchronize the satellite's spin to its orbital 
period and circularize orbits) at larger semi-major axes. 
The exploration of evolutionary processes that can ex- 
plain these spin-orbital characteristics will be discussed 
in the bulk of this paper. 

Plausible 1—a uncertainties for semi-major axes and 
eccentricities are compiled from published values (when 



available) and listed in Table [T] along with the adopted 
values. When published uncertainties are not available, 
we obtain uncert ainties by determining their 1—cr confi- 
dence levels (e.g. lCash|[T976tlPress et al.lll992f ). For each 
considered parameter, we hold it fixed at a range of val- 
ues while simultaneously fitting for all other parameters. 
Since only one parameter is held fixed at a time, a I— a 
confidence region is prescribed by the range of solutions 
which exhibit chi-square values within 1.0 of the best-fit 
chi-square. 

All of these well-characterized NEA systems are ob- 
served by radar and have varying degrees of observa- 
tional coverage and quality. We have high confidence 
in the datasets for binaries 1999 KW4, 2000 DP107, and 
1991 VH as weU as triples 2001 SN263 and 1994 CC. 
These systems have extensive, high signal-to-noise mea- 
surements on ^10 epochs over ^2 weeks, or have been 
observed on 2 separate apparitions. We have moderate 
confidence for 2002 CE26, Didymos, 2004 DC, and 2003 
YTl, which all have ~50— 150 measurements over 4—14 
days. These high and moderate confidence datasets com- 
prise the well-characterized sample of NEA binaries and 
triples listed in Table [T] We also mention another radar- 
observed binary designated 1998 ST27 with a n asyn- 
chronously rotating satellite (iBenner et al.l 120031 ) due to 
its uniquely large separation. 1998 ST27 has a low con- 
fidence dataset with fewer than 40 measurements and 
inconsistencies in orbit solutions, and thus we can only 
determine a lower bound on its semi-major axis of >12 
primary radii or >5 km. Its actual semi-major axis may 
be much higher. Its large separation makes 1998 ST27 
the widest NEA binary discovered so far. 

Since all of these NEA systems in Table [1] are char- 
acterized by radar data, this sample is biased towards 
binaries and triples discovered through radar techniques, 
which typically require close approaches with Earth of < 
0.1 AU. Close approach data, including the most recent 
approach near the time of radar observations and the cur- 
rent minimum orbital intersection distance (MOID) with 
Earth, arc hsted in Table [H The MOID describes the 
minimum distance between two elliptical orbits and dis- 
regar ds the positions of the bodies in their orbits (jSitarskil 
|1968[ ). and is valid as long as the osculating orbital el- 
ements approximate the actual orbits. These orbital el- 
ements certainly change over long periods of time and 
therefore close approaches to Earth less than the current 
MOID could have occurred in the past. 

1.2. Dynamical and Collisional Lifetimes 

NEAs are short-lived and have d ynamical lifetimes o n 
the order of a few million years (jBottke et al.l [200l . 
These average lifetimes represent how long they can sur- 
vive in near-Earth space before plunging into the Sun, 
getting ejected from the Solar System, or colliding into a 
planet. Due to these short NEA lifetimes, the near-Earth 
population is continually replenished by small (^10 km) 
MBAs. These small MBAs migrate into near-Earth space 
through unstable main belt regions that are permeated 
by strong resonances with Jupiter and Saturn. Reso- 
nance capture is enabled by radiation forces and colli- 
sions, which bring MBAs into these unstable main belt 
regions. While in the main belt, these asteroids have col- 
lisional lifetimes that are dependent on their sizes. For 
example, an asteroid 500 meters in radius will have a col- 
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Figure 1. This plot shows the well-characterized NEA 
systems in eccentricity and semi-major axis (in primary 
radius) space, whose orbital parameters are taken from 
Table [TJ Asynchronous satellites are represented by un- 
filled circles and synchronous satellites are marked with 
filled circles. 



hsional hfetime of -3.8 x 10^ years (jBottke et al.l l2005D. 

Clearly, collisional lifetimes (while in the main belt) 
are much greater than dynamical lifetimes (while in near- 
Earth space). It is quite plausible that many observed 
NEA systems had their satellites formed while still in 
the main belt, and so their total lifetimes as a binary or 
triple will be dominated by their prior collisional lifetime. 
An important implication of binary /triple formation in 
the main belt is that some evolutionary processes have a 
longer period of time within which they may occur. The 
examination of these evolutionary processes constitute 
the bulk of this paper. 

1.3. Binary and Triple Formation 

Binary and triple asteroids form t hrough the generally 
accepted rot ational fiss ion model (IMargot et all 120021 : 
iPravec et all [2006; Walsh et al.ll2008| ). The likeh^spin- 
up m echanism is the thermal YORP effect ([Rubincaml 
[2000h . which causes mass shedding from the primary. In 
this subs ection, we introduce t he po st-fission dynamics 
model of iJacobson fc Scheered (|2011| ). and in particular 
we evaluate its relevance in explaining any of the ob- 
served spin-orbital characteristics of NEA binaries and 
tr iples. 

IJacobson fc ScheeresI ()201in describe the immediate 
(<1000 years) dynamics following rotational fission, 
which include YORP (to spin up the initial body, and 
this is the only time a non-gravitational force is incor- 
porated in their model), secondary fission (a satellite is 
rotationally accelerated and then fissions to create an- 
other satellite), tri-axial gravitational potentials, tides, 
and solar gravitational perturbations. The immediate 
result after initial fissioning of the primary is a chaotic bi- 
nary, and subsequent evolutionary processes are mainly 
determined by mass ratio (mass of secondary/mass of 
primary) . 

Binaries with high mass ratios (> 0.2) do not undergo 
secondary fission and instead experience tidal dissipa- 
tion to become doubly synchronous. This may lead to 
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a contact binary state if BYORP contracts the orbit. 
Low mass ratio binaries will mostly disrupt unless they 
are allowed to experience secondary fission, which cre- 
ates an initially chaotic triple system. The chaotic triple 
can stabilize and become a binary by ejection or colli- 
sion of a satellite, which can lower the system's angular 
momentum and energy. Resulting binaries that are sta- 
ble after 1000 years of evol ution are shown in Figure [2] 
(jJacobson fc Scheeresll2011t ). 

In Figure [21 we point out two scenarios of inter- 
est. First, early tidal evolution can lead to a com- 
monly observed type of binary: a synchronous satellite 
with a separation of ~4 primary radii and an eccentric- 
ity of ^0.3—0.4 that will continue to damp by tides. 
Second, some satellites have large separations (up to 
'--^lO— 16 primary radii) from their primary with corre- 
spondingly large eccentricities (up to ^0.6—0.8). These 
wide primary-satellite pairs may explain observed satel- 
lites located at large semi-major axes, although the ec- 
centricities predicted by the post-fission dynamics model 
are larger than any observed value. 

Triple formation can occur through secondary fission 
or another round o f YOR P-induced primary fission. 
iJacobson &: ScheeresI ()2011[ ) include secondary fission in 
this model and their simulations produce no stable triple 
systems after 1000 years of evolution. They do not 
model additional rounds of YORP-induced primary fis- 
sion, which remains a plausible explanation for the for- 
mation of triple systems. We add the hypothesis that 
triples form by first creating a wide binary (such as those 
seen in post-fission simulations) through primary fission, 
and then a closer-in satellite is formed during a subse- 
quent round of primary fission. This possibility is also 
supported by observations of wide binary 1998 ST27 that 
has a separation (>12 primary radii) consistent with the 
outer satellite's separation (^13 primary radii) in triple 
2001 SN263. 

Other studies discussing spin-up fissioning either do 
not mention or only provide scant information about 
the resulting spin-orbital parameters of a newly-formed 
satell ite (e.g. lScheereil2007tlHolsapplell2010l:[Walsh et al.l 
[20081 ). The fissioning model of iWalsh et al.l (|2008l ) does 
not attempt to simulate post-fission dynamics and so 
does not include tidal interactions. In their simulations, 
when satellites grow to 0.3 primary radii, their separa- 
tions are 2—4.5 primary radii and their eccentricities are 
<0.15. This range of eccentr i cities is lower than that 
found by IJacobson fc ScheeresI (|2011[ ) and shown in Fig- 
ure [H 

In essence, the immediate (<1000 years) dynamics en- 
suing from this formation scenario provide a pathway 
for the creation of wide binaries such as 1998 ST27 and 
potentially the outer satellites in triple systems. How- 
ever, the eccentricities from post-fission dynamics are too 
high compared to observed eccentricities, and some of the 
simulated binaries have small semi-major axes of ~2 pri- 
mary radii that are smaller than any observed separation. 
Moreover, the spin states of just-formed satellites will be 
asynchronous. Clearly, if multiple systems are formed 
by rotationa l fission and follow th e post -fission dynam- 
ics model of IJacobson fc ScheeresI (|201lD . there will be 
additional processes that evolve the systems, and the ex- 
ploration of these processes constitute the bulk and re- 



mainder of this paper. 
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Figure 2. Results f rom post-fission simulations by 
IJacobson fc ScheeresI ()2011[ ) are shown after 1000 years 
of evolution (Seth Jacobson, personal communication, 
2011). 

2. TIDAL EVOLUTION 

In this section, we investigate evolutionary processes 
due to tidal effects, and whether tides can explain the 
observed spin-orbital characteristics of satellites in NEA 
systems. The fastest tidal evolution process (in the ab- 
sence of other spin-modifying forces) is the synchroniza- 
tion of a satellite's spin to its orbital period due to the 
tides raised on the secondary by the primary ([Goldreichl 
119631 iGoldreich fc Saril |2009| ). After synchronization, 
tidal evolution continues by modifying the semi-major 
axis and eccentricity of the mutual orbit. It is a compet- 
ing process between the opposing effects of tides raised on 
a primary (which increase both the semi-major axis and 
eccentricity) and those raised on the secondary (which 
cause negligible changes in the semi-major axis, but de- 
crease eccentricity). Tidal evolution ends after the orbit 
has circularized and the primary's spin has synchronized 
to the mutual orbit period, resulting in a doubly syn- 
chronous system. 

Models of tidal evolution for asteroids are governed 
by two dimensionless parameters: the effective rigid- 
ity jl and the tidal dissipation factor Q. The non- 
dimensional effective rigidity /i is a function of the 
body's internal properties such as density p and ra- 
dius R, and is defined as follows for a monolith: jl = 
{l9fj,)/(2gpR) where g represents self-gravity and fi is 
the b ody's rigidity or shear mo dulus with units of pres- 
sure (jMurrav fc Dermottlll999l ). The effective rigidity is 
related to the tidal Love number fc, where k — 1.5/(l+/i) 
for a homogeneous solid body. The tidal dissipation fac- 
tor Q is a quality factor defined as Q = {2ttEo)/{AE), 
which d escribes the body's effectiveness at dissipating 
energy (jMurrav fc DermottI [19991 ). If we consider the 
body's response to tidal oscillations as a harmonic oscil- 
lator, Eq represents the peak energy stored during a cycle 
and AE is the energy dissipated over a cycle. Unfortu- 
nately, for small asteroids there are substantial uncer- 
tainties in crucial parameters jl (or k) and Q, in addition 
to other poorly known effects, such as the frequency de- 
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pendence of these two quantities, their apphcabihty to 
porous bodies, or even our abihty to capture tidal pro- 
cesses with two ideahzed numbers. 

2.1. Eccentricity Evolution 

Post-fission eccentricities (^Section ll.3p are significantly 
higher than observed eccentricities for some NEA sys- 
tems, and so here we explore the effects of tides in mod- 
ifying eccentricity. Several models have been considered 
to facilitate calculations of tidal evolution in asteroids. 

In the monolith model, asteroids are idealized as uni- 
form bodies with no voids. In this idealization, the effec- 
tive rigidity /imono is inversely proportional to the square 
of the asteroid's radius R: /jmn nn oc or k 



For all models, the eccentricit y evoluti o n for 
a synchronous sate llite's orbit ( GoldreichI 119631 : 
IGoldreich fc^ari|[2009l) is 



To arrive at numerical values, IGoldreich &: Saril ()2009l ) 
used the Moon's radi us and Love number k of ^0.025 
([Williams et al.l 120081 ) and obtained /Imono ^ 2 x 10^ 
(1 km/i?)^. The monolith tidal model has been used 
to estimate the relative strengths of the components in 
the Kalliope-Linus binary system. Since its mutual orbit 
is found to be near-circular, comparison of the relative 
rates of eccentricity excitation and damping constrains 
the relative values of ilQ for the primary and secondary: 
iMargot fc Brownl ()2003l ) found that jlQ for the secondary 
is smaller than that of the primary. 

In the rubble pile model, asteroids arc idealized as grav- 
itational aggregates, i.e. composed of smaller elements 
held together by gravity only. This assumption is moti- 
vated by the low observed d ensities of NEA sys tems and 
the observed "spin barrier" (jPravec et al.ll2002l ). We wiU 
c onsider two separate rub ble pile models. 

IGoldreich fcSadI (|2009h propose that the relationship 
between a rubble pile's effective rigidity /Irubbio and a 
monolith's effective rigidity /imono of comparable com- 
position an d size is simply Mrnhhip ^ lOVAmono- Thus, 
this model (jGoldreich fc Sarl2009[ ) defines a rubble pile's 
effective rigidity as jlmhUe ~ lO^V^ (1 km/i?) or Love 
number A:iubbic ~ 10~^ km). If we assume a com- 

mon density and tidal quality Q factor between the pri- 
mary and secondary, the rubble pile model adopted by 
IGoldreich fc Sari (2009) gives the ratio of the rates of 
eccentricity excitation to damping as 19/28 (irrespective 
of component sizes) for a system with a synchronized 
secondary. Therefore, in this model the eccentricity will 
likely damp for such systems. In cases where there are 
density or Q differences between the components, it is 
possible t he eccentricity ca n grow . 

iJacobson fc ScheeresI (j2011bl ) find that the 
IGoldreich fc Saril (120091 ) model agrees reasonably 
well for data of binaries with primary radii of ^2 km. 
but not for systems with very small primary radii. 
This discrepancy can be resolved using a different Love 
number dependence on size, and IJacobson fc ScheeresI 
(|2011b[) describe a power law fit to the log arithmic data 
of kno wn synchronous binaries compiled bv lPravec et al.l 
(|2006[) . They find that /irubbic - 6 x 10" {R/1 km) 
or fcrubbiG ~ 2.5 X 10~^ (1 km/R). Thus, if we apply 
this model for binary components of common density 
and tidal quality factor Q, we find that the ratio 
of the rates of eccentricity excitation to damping is 
(19iZs^)/(28i?p^), where Rp is the radius of the primary 
and Rs is the radius of the secondary. Thus, eccentricity 
will also likely damp in this model. 



de 57 kp Ms f Rp 



dt 



]pMp 



21 ks Mp (Rs . 

Yq-sWs [- ) W 



where there are two competing terms corresponding to 
eccentricity excitation and damping, respectively. This 
equation is a function of tidal Love number k, tidal qual- 
ity factor Q, mass M, radius R, semi-major axis a, mean 
motion n, and eccentricity e. The subscripts p and s 
represent the primary and secondary, respectively. 

Using Equation [1] we calculate the circularization 
timescales rjamp of all orbits in our NEA sample given 
their current orbital and physical parameters (Table [ij . 
These timescales are presented in Table [2] In our cal- 
culations, we use an assumption of Q 10—10^, which 
is reasonable for small rocky bodies. We calculate rjamp 
ac cording to both rubble p ile models: the k oc R model 
of IGoldreich fc Saril (|2009|) and the k (x R-^ model 
of Ijacobson fc Scheered (1201 Ibf ). With these assump- 
tions, we obtain Tdamp by calculating the 1/e damping 
timcscalc (here e is Euler's constant) for eccentricity by 
numerically integrating Equation [1] 

For the cases of 1999 KW4, 2001 SN263 #1, 
1994 CC #1, and 1994 CC #2 the tidal Love number 
model of lGoldreich fc Saril (120091) gives a larger de/dt for 
the excitation term than for the damping term (which oc- 
curs for these cases due to density differences between the 
primary and satellite), and so their eccentricities would 
theoretically be predicted to increase. However, most of 
these satellites are observed to have circular orbits, which 
indicates that either the densities are incorrect, there are 
Q differences between the primary and satellite, and/or 
the tidal model is incorrect. For these cases, we suspect 
that the satellite densities are in error. If we assume that 
all components have a density equal to that of the pri- 
mary when we compute de/dt, then we find that their 
eccentricities will damp as predicted by theory. 

The eccentricity damping timescales calculated from 
Equation [T] and presented in Table [2] show that there 
are differences of several orders of magnitude between 
times cales calculated from th e tw o different tidal mod- 
els oflGoldreich fc Saril ()2009[ ) and IJacobson fc ScheeresI 
(|2011b( ). As expected, satellites with larger semi-major 
axes such as the outer satellites in 2001 SN263 and 
1994 CC will have mutual orbits that take longer to cir- 
cularize, and closer-in systems will have shorter damping 
timescales. The total possible lifetimes of each system, 
which includes collisional and dynamical lifetimes, are 
also listed in Table O For 1994 CC's outer satellite, its 
orbit has a damping timescale that is clearly greater than 
its total possible lifetime. If this satellite formed from 
a po st-fission dynamics scenario (jJacobson fc ScheeresI 
120 lit) with a high post- fission eccentricity (~0.8), the sys- 
tem cannot evolve to the currently observed eccentricity 
of 0.19 by tides within its lifetime; this remains true even 
if there are Q differences between the primary and sec- 
ondary (i.e. let Qp = 100 and Qs = 10) or if the satellite's 
size is underestimated (its r adius can be as large as 55 
meters: iBrozovic et all 120111 ). If 1994 CC's outer satel- 
lite is to be explained by post-fission dynamics and tides, 
then either fissioning can deliver a far-flung satellite with 
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lower eccentricities or the equations that ideahze tidal in- 
teractions are not sufficient. For 1999 KW4, 2002 CE26, 
Didymos, 1991 VH, and 2001 SN263 #1, their ranges 
of damping timescales are less than their total possible 
lifetimes and this suggests that they should be close to 
circular, which is corroborated by the fact that their ob- 
served e are all less than 0.1. For most NEAs, we cannot 
draw firm conclusions since these calculations are very 
dependent on crucial parameters such as k and Q, which 
are poorly constrained for small asteroids; their damping 
timescales can be either greater or less than their total 
lifetimes. 

To summarize, in this subsection we have considered 
the effect of tides in modifying eccentricity. Current 
models show that tides will damp eccentricity, but in- 
creasing eccentricity is possible if the primary and sec- 
ondary have different densit y and/or Q values. The post- 
fission dynamics model of iJacobson fc ScheeresI (|2011[ ) 
tends to produce high-e in distant satellites, and we find 
that for all satellites except the outer satellite of 1994 CC, 
it is possible that tides can damp their post-fission ec- 
centricities to observed eccentricities within a collisional 
lifetime. 

2.2. Satellite's Spin Evolution 

We consider a satellite's spin evolution due to tidal 
torques (on the tide-generated bulge as well as on the 
permanent deformation) and the radiative YORP effect. 
Spin evolution is important in our investigation on the 
origin of the observed spin states of NEA satellites as 
well as its consequences for spin-orbit synchronization- 
dependent processes such as BYORP. 

In the absence of other perturbations, spin-orbit syn- 
chronization of the sat ellite is the fastest tidal process 
(jGoldreich fc Sarill2009t) according to 



— = 5.sgn(n-a..)-Gp,(^-j 



(2) 



where lo is the spin rate, n is the mean motion, k is 
the Love number, Q is the tidal dissipation factor, G is 
the gravitational constant, p is the density, R is the ra- 
dius, and a is the semi-major axis. Subscripts p and s 
denote the primary and satellite, respectively. All ob- 
served asynchronous satellites in our sample are fast ro- 
tators that spin faster than their orbital motions and so 
diOg/dt will be negative. We calculate the synchroniza- 
tion timescales Tgync to despin from a breakup period 
(^--^2.33 hours) to its orbital period according to Equa- 
tion [21 and present them in Table [31 These despinning 
timescales are shorter than the total lifetimes for most 
satellites and are shorter than the eccentricity damping 
timescales for all satellites (Table [5]). These synchroniza- 
tion timescales will likely be affected by YORP, which 
may help speed up synchronization or slow it down. 

Accordingly, we investigate the effect of YORP during 
a satellite's spin evolution. For both cases of an initially 
asynchronous or synchronous satellite, we ask, how does 
a satellite's rotation evolve when the semi-major axis 
changes? We do not assume that the orbit migration 
is dominated by any particular mechanism. Perturba- 
tions such as tides or BYORP can increase a satellite's 
semi-major axis, and BYORP can also decrease the semi- 
major axis. 



2.2.1. Initially Asynchronous Satellite 

First, we examine an initially asynchronous satellite, 
perhaps newly formed. All satellites will go through this 
stage. Its spin evolution will be affected by YORP as 
well as the torque due to the satellite's tidal bulge. The 
restoring torque due to the satellite's permanent defor- 
mation is not applicable in this case because this torque 
averag es out unless the satellite is in a spin-orbit reso- 
nance (|Murrav fc Dermott. 19991 

We apply the torque equations due to tides raised 
on the satellite by t he primary and YORP using 
the formalism gi v en in iMurrav fc DermottI (I1999D and 
ISteinberg fc Sari ([20lTh for the magnitudes of the 
torques N: 

3 fc, GM^Rl 

Midn = ^ ^ 



2Qs 



YORP 



6cd| 



(3) 



(4) 



1 



where G, c, and Lq are the gravitational constant, the 
speed of light, and the solar luminosity, respectively. The 
heliocentric parameters include semi-major axis cIq and 
eccentricity e©. Tidal parameters include the Love num- 
ber ks and dissipation factor Qs, where the subscript s 
denotes that these quantities are for the satellite. Mp is 
the primary's mass, a is the se mi- major axis, and Rs i s 
the sateUite's radius. Following [Goldreich fc Saril (|2009[) . 
we include a YORP reduction factor /y that can be pos- 
itive or negative, and accounts for a reduction from its 
maximum possible value. This factor is necessary since 
the incoming radiation is not completely absorbed and 
reemitted tangentially along the satellite's equator. 

The tidal torque will try to establish synchronization 
between the secondary's spin and its orbital period, and 
the YORP torque will perturb the secondary's spin in 
cither direction. Their relative contributions are depen- 
dent on the semi-major axis. If the torques act in the 
same direction, the satellite's spin will synchronize with 
its orbital period. If the torques are competing, spin syn- 
chronization is not guaranteed. The critical semi-major 
axis ac,tide at which the magnitudes of these torques are 
equal is 



ac,tido = (3i?sMpdo) 



1/3 



KGcJl 



1/6 



Qs\fY\L, 



(5) 



where the constants and variables are the same as defined 
for Equations [3] and m 

Satellites with orbital distances less than this critical 
semi-major axis Oc^tidc will become synchronous, whether 
or not YORP and tides act in the same direction. For or- 
bital distances greater than the critical semi-major axis, 
satellites will likely remain asynchronous if tides and 
YORP act in opposite directions. We calculate a^tido for 
all systems in our sample (Table[T|), using Qs = 10— 100 
and /y ^ 5xl0~^as suggested by observations of 
(54509) YORP (fo rmerly 2000 PH5) ()Lowrv et al.|[2007l: 
iTavlor et all 120071 ). A list of ac,tido is given in Table H 
Uncertainties in the tidal Love number fc, and dissipa- 
tion value Qs produce a range of critical semi-major axes 
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Table 2 

Tidal Timescales for Eccentricity Damping 



System 


a/Rp 


e 


''"lifetime 


(yr) 


(fc oc 


'''damp 


(yr) 


(fc oc 


R-^ 


t 


2000 DP107 


6.6 


0.01 


3. 


,36 


X 


10« 


1 X 10* 


- 1 X 10'-* 


4 


X 


10^ 


- 4 


X 


10® 


1999 KW4 


3.9 


0.008 


4, 


,32 


X 


10* 


3 X 10^^ - 


- 3 X lO'^t 


3 


X 


10"* 


- 3 


X 


10^ 


2002 CE26 


2.8 


0.025 


7. 


,03 


X 


10* 


1 X 10" 


- 1 X 10** 


4 


X 


10-* 


- 4 


X 


10^ 


*2004 DC 


4.4 


0.30 


2. 


19 


X 


10* 


2 X 10** 


- 2 X 10^ 


3 


X 


10-* 


- 3 


X 


10^ 


*2003 YTl 


7.1 


0.18 


3. 


,94 


X 


10** 


1 X 10^ - 


- 1 X 10^° 


1 


X 


lO'' 


- 1 


X 


10^ 


Didymos 


3.0 


0.04 


3. 


,36 


X 


10* 


5 X 10® 


- 5 X lO'^ 


4 


X 


10^ 


- 4 


X 


10* 


*1991 VH 


5.4 


0.06 


4. 


12 


X 


10« 


2 X 10^ 


- 2 X 10** 


2 


X 


10^ 


- 2 


X 


10® 


2001 SN263 #1 


2.9 


0.016 


6. 


06 


X 


10* 


3 X lO'^ - 


- 3 X lO'^t 


4 


X 


10* 


- 4 


X 


105 


*2001 SN263 #2 


13 


0.015 


6. 


,06 


X 


10* 


4 X 10** 


- 4 X 10^ 


4 


X 


10^ 


- 4 


X 


10* 


1994 CC #1 


5.6 


0.002 


2. 


,96 


X 


108 


5 X 10** - 


- 5 X lO^t 


1 


X 


10® 


- 1 


X 


10^ 


*1994 CC #2 


20 


0.19 


2. 


96 


X 


10« 


5 X 10^2 - 


5 X W^^t 


2 X 10** - 


- 2 X lO^" 



For each system's satellite, we list the adopted values for the observed semi-major axis 
a/Rp and eccentricity e, its t otal possible lifetime "nifctimc (coUisional lifetime plus dynamical 
lifetime) (jBottke et al.ll2"005l ). and the eccentricity damping timescale rjamp due to tides for 
synchronous satellites. We include two tidal models with a different size dependence of the 
Love number k. For each model, there is a range of values due to adopted values for the tidal 
dissipation factor Q ~ 10—10^. 
* Known asynchronous satellites 

X Systems in which we assumed that all components had density equivalent to that of primary 
IGoldreich fc Saril ||2009|) 
' iJacobson fc ScheeresI I l2011bl ') 



Table 3 



Tidal Timescales and Distances for Spin Synchronization 


System 


a/Rp 


e 


'''sync 


(yr): 




1c,tido 


/Rp: 




Rh (km) 


Clc.porm (km) 








(k oc R)^ 


(fc oc 




(fc OC R]^ 


{k oc 








2000 DP107 


6.55 


0.01 


7 X 10® - 7 X 10^ 


6 X 10"* 


- 6 X 10^5 


2.88 - 4.23 


6.33 


- 9.28 


86 


1091 


1999 KW4 


3.86 


0.008 


2 X 10^ - 2 X 10® 


4 X 10=* 


- 4 X 10"' 


2.80 - 4.12 


5.36 


- 7.86 


70 


1703 


2002 CE26 


2.78 


0.025 


6 X 10"* - 6 X lO'"^ 


6 X 10^ 


- 6 X 10=* 


2.80 - 4.11 


6.14 


- 9.01 


514 


4584 


*2004 DC 


4.41 


0.30 


3 X 10® - 3 X lO'^ 


1 X 10=* 


- 1 X 10"* 


1.41 - 2.06 


5.24 


- 7.65 


44 


189 


*2003 YTl 


7.15 


0.18 


1 X lO'^ - 1 X 10** 


4 X lO"' 


- 4 X 10^ 


2.67 - 3.93 


6.62 


- 9.71 


113 


1488 


Didymos 


2.95 


0.04 


8 X 10"* - 8 X lO'"^ 


2 X 10^ 


- 2 X 10=* 


2.30 - 3.38 


6.35 


- 9.33 


109 


891 


*1991 VH 


5.43 


0.06 


1 X 10® - 1 X lO'^ 


3 X 10'* 


- 3 X 10=^ 


3.40 - 5.00 


6.38 


- 9.37 


105 


1987 


2001 SN263 #1 


2.92 


0.016 


1 X lO"^ - 1 X 10® 


2 X 10=* 


- 2 X 10-* 


3.40 - 4.99 


6.46 


- 9.48 


343 


5415 


*2001 SN263 #2 


12.8 


0.015 


7 X lO'^ - 7 X 10** 


8 X 10® 


- 8 X lO'^ 


5.16 - 7.58 


7.43 


- 10.9 


343 


5532 


1994 CC #1 


5.58 


0.002 


2 X lO'^ - 2 X 10** 


3 X 10'' 


- 3 X 10^ 


2.03 - 3.00 


6.16 


- 9.06 


86 


911 


*1994 CC #2 


19.8 


0.19 


3 X 10^'' - 3 X 10" 


2 X lO'^ 


- 2 X 10** 


1.71 - 2.52 


5.81 


- 8.55 


86 


553 



For each system's satellite, wc list the adopted values for the observed semi-major axis a/Rp and eccentricity e, the tidal spin 
synchronization timescale Tgync starting from the breakup rate in the absence of other effects such as YORP, the critical semi-major 
axis fic,tido/^p within which an initially asynchronous satellite can achieve synchronization against the effects of YORP torques that 
compete with tidal torques. Hill radii Rh, and an upper limit to the critical semi-major axis ac,perm at which an initially synchronous 
satellite would break spin lock due to the effects of YORP. The critical semi-major axes are calculated using fy ~ 5 X 10~*. When 
applicable, we include 2 tidal models with a different size dependence of the Love number k. For each model, there is a range of values 
due to adopted values for the tidal dissipation factor Q ~ 10—10^. 
* Known asynchronous satellites 

IGoldreich & Sari (JO O S) 

' iJacobson fc Scheerej I l2011bl 1 
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per NEA binary in our sample, and different values for 
fy would result in an even broader range. 

We now discuss the agreement between our calcula- 
tions of article and the observed spin states (given in Ta- 
ble of NEA satellites in our sample. Two satellites, 
2002 CE26's secondary and 2001 SN263's inner satel- 
lite, have current semi-major axes less than the range of 
possible Octide values calculated in this table, and they 
are both observed to be synchronous. Two satellites, the 
outer satellites of 2001 SN263 and 1994 CC, have current 
semi-major axes larger than their range of possible ac,tido 
values, and they are both observed to be asynchronous; 
this suggests that tides and YORP act in opposing direc- 
tions for these two satellites. All other satellites whose 
current semi-major axes may or may not be larger than 
their Oc.tido include synchronous and asynchronous ro- 
tators. The dominance of YORP for satellites such as 
2001 SN263 #2 explains how it could have a tidally- 
circularized orbit yet be asynchronous. 

Further observations of asyn chronous satellites can 
help constrain which tidal m odel (jGoldreich &: Sarl2009l : 
iJacobson fc Scheeresll20l"Tbl ) best captures the behavior 
of small bodies, assuming that tides and YORP are the 
dominant processes affecting satellite spin states. For 
this test we focus on asynchronous satellites because in 
that case the balance of tidal and YORP torques requires 
their observed semi-major axes to be larger than their 
computed ac.tidc values. Comparison of the observed a to 
the predicted flctido for a number of asynchronous satel- 
lites of different sizes may reveal which tidal model is 
more appropriate. We encourage observations of asyn- 
chronous satellites to enable this test. Unfortunately the 
same test cannot be applied to the more numerous syn- 
chronous satellites because their spin can be explained 
by tides and YORP torques acting in the same direction 
with the observed a < a^tide or a > Oc^tide, or by tides 
and YORP acting in opposite directions with a < Octido- 

2.2.2. Initially Synchronous Satellite 

We now consider the case of an initially synchronous 
satellite. Its spin evolution will be affected by YORP, 
the satellite's tidal bulge, and additionally, a restoring 
torque due to the satellite's permanent deformation. The 
permanent quadrupole moment of the satellite plays a 
role for spin states in spin-orbit resonance and is thus 
applicable when considering the spin evolution of a syn- 
chronous satellite. The expression for the magnitude of 
the torque due to perma nent deformation is given by 
IMurrav fc DermottI (jl999D as 



(B_A)^sin(2V') 



(6) 



where B and A are the satellite's equatorial principal 
moments of inertia, G is the gravitational constant. Alp 
is the mass of the primary, a is the semi-major axis, and 
ip is the amplitude of the libration. 

We compare the tidal torque iVtido to the permanent 
deformation torque A^pcrm for the case of 1999 KW4, 
which is currently the only NEA binary with published 
shape information of the secondary. Using 1999 KW4's 
(B — A) value for its secondary, ■0 ~ 7r/4 (maximum am- 
plitude), values for Q from 1 to 100, and b oth tidal mod- 
els (|Goldreich fc Sari,2009; Jacobson fc Sclieeres„2011bi ). 



we find that A'porm dominates over iVtide by at least six 
orders of magnitude. Therefore, for the case of an ini- 
tially synchronous satellite that we examine here, we will 
consider only the permanent deformation torque (and 
not the torque on the tidal bulge) as well as the YORP 
torque. The permanent deformation torque will seek 
to maintain spin-orbit synchronization, and the YORP 
torque will attempt to spin the satellite in either direction 
away from synchronization. Balance of the permanent 
deformation and YORP torques yields a critical semi- 
major axis ac,porm given as 

1/3 



I 9{B- A)GMp sin(2?/')cd|, /I 



-''C,perm 



(7) 

where the constants and variables are the same as defined 
for Equations m and ini 

In the absence of significant eccentricities, this critical 
semi-major axis ac,pcim separates regions where synchro- 
nization can be maintained and where synchronization 
can be broken. We seek an upper limit for Ocperm to 
determine if synchronization can be maintained to dis- 
tances as large as the Hill radius, which is a requirement 
for synchronization-dependent processes such as BYORP 
to create asteroid pairs. In Table [3l we present calcula- 
tions of an upper limit ac,perm for all satellites in our sam- 
ple. We adopt a maximum libration amplitude ip ~ 7r/4 
and calculate {B — A) for each satellite by scaling from 
the {B — A) of one of the mos t elongated NEAs known, 
Geographos (jOstro et al.lll995D : 



[B - A)^ 



(8) 



where subscripts x and Geo represent the considered 
satellite and Geographos, respectively. The mass is 
M and the equivalent radius (if it were spherical) is 
R. The equivalent rad ius of Geographos is ^^1.28 km 
(|Hudson fc Ostrol[T99l . and if we assume a typical rub- 
ble pile density of 2 g cm"'^, its mass is ~^1.8 x lO^'^ 
kg. Using a triaxial ellipsoid assumption to calculate its 
moments of inertia, we find {B — A)Geo/(A^Goo-RGoo) to 
be ^0.56. Similarly as we did earlier for calculations of 
tie, tide, we adopt fy ^ ^'x 10"''. These values are used in 
the calculation of an upper limit Oc porm shown in Table 

El 

Tablc[3]shows that for all NEA systems, ac,perm is much 
larger than the Hill radius. This suggests that once a 
satellite has synchronized its spin to its orbital period, it 
can remain so out to far distances even in the presence 
of YORP. Therefore, all observed asynchronous satellites 
(2004 DC, 2003 YTl, 1991 VH, and the outer satellites 
of triples 2001 SN263 and 1994 CC) have probably never 
been synchronous unless they experienced a planetary 
encounter that disrupted their synchronous spins. 

An alternative proposal w as recently put forth by 
IJacobson fc ScheeresI (|2011a[ ). who suggest that longi- 
tude librations of the satellite will occur about a direction 
that is not aligned with the line connecting the central 
body and the satellite. They hypothesize that this an- 
gular offset becomes increasingly significant as the orbit 
expands (due to tides and BYORP), and eventually re- 
sults in breaking a satellite's spin lock. If this model is 
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correct, then synchronization cannot be maintained to 
far distances. However this model seems to produce sub- 
stantial angular offsets only for satellites with very small 
moment of inertia differences ((B — A)/C), and there- 
fore docs not appear to be effective for the overwhelming 
majority of satellites. 

To summarize, we find that a newly post-fissioned, 
asynchronous satellite can become synchronous within 
its lifetime and this can explain all observed synchronous 
and circular binaries and inner satellites of triples. Asyn- 
chronous outer satellites in both triples are unable to syn- 
chronize by tidal torques because of the larger influence 
of YORP, and this may also be the explanation for asyn- 
chronous binaries 2004 DC, 2003 YTl, and 1991 VH. 
The dominance of YORP over tides at large distances 
may explain why the outer satellite in 2001 SN263 re- 
mains asynchronous despite having an orbit that appears 
to have tidally circularized. 

2.3. Semi-major Axis Evolution 

The post-fissio n dynamics model of 

iJacobson fc ScheeresI ()2011f ) shows how satellites 
can be deposi ted a t large separations up to ^16 primary 
radii (Section II. 3p and can explain the outer satellite in 
2001 SN263 (at ~13 primary radii) and perhaps also 
the outer satellite in 1994 CC (at ^20 primary radii). 
Here we investigate if another mechanism, namely tidal 
evolution, can bring satellites from closer-in to wide 
orbits. 

To explore this possibility, we study three cases of 
wide orbits: 2001 SN263 #2, 1994 CC #2, and we 
also consider a hypothetical wide binary modeled af- 
ter 1998 ST27, which has a primary radius of 0.42 km 
and a lower limit on its observed semi-major axis of 
5 km. For these three systems, we calculate the tidal 
timcscales for the semi-major axis to increase from one 
pri mary radius to its observed value. For both tidal mod- 
els (IGoldreich fc Saril[2009l : IJacobson fc ScheeresI [20TTbh 
and using Q = 10—100, we find that the tidal 
timescales for such increases in semi-major axis are 
-10^-101° years for 2001 SN263 #2, ~10i°-10i3 years 
for 1994 CC #2, and -10^-10" for the hypothetical 
wide binary modeled after 1998 ST27. The total possi- 
ble lifetimes of these satellites (see Table [21 for the hy- 
pothetical binary, its lifetime is ~3.4 x 10^ years) are 
shorter than their tidal expansion timescales. This sug- 
gests that tides cannot account for the wide orbits of 
these satellites, if the tidal model is correct. Other influ- 
encing factors such as "tidal saltation," the idea of mass 
lofting from the primary and briefly entering into orbit to 
transfer orbital angular momentum to the satellite before 
falling back dow n to the primary's surface, may speed up 
tidal expansion (jHarris et al.l[200l . 

Next, wc investigate if tides can explain the wide orbits 
by using different assumptions of interior properties. As- 
suming tidal dissipation Q values of 10— 100, we calculate 
their material properties in order for their semi-major 
axes to increase from one primary radius to the currently 
observed semi-major axis within their NEA-specific life- 
times. We find that in order for tides to be responsible 
for the increase in semi-major axis within the considered 
system's lifetime, we will need to invoke the following 
values for kp-. -0.00003-0.00030 for 2001 SN263 #2, 
~0.0054-0.0545 for 1994 CC #2, and ~0.00028-0.00285 



for the hypothetical binary modeled after 1998 ST27. 
The required k value for 1994 CC #2 is prohibitive and 
is almost as larg e as or larger than the Moon's k of 0.025 
(I Williams et aLl[2008h . 

In summary, tidal evolution cannot explain the semi- 
major axes of widely-separated systems without invok- 
ing unusual material properties, lower Q values than as- 
sumed, a different tidal model, or a combination of these 
factors. Our analysis here supports the idea that post- 
fission dynamics may be largely responsible for some of 
the wide orbits observed in NBA systems. 

3. BYORP EVOLUTION 

The BYORP effect occurs for a synchronously rotat- 
ing satellite with an asymmetrical shape. A synchronous 
satellite has permanent leading and trailing hemispheres, 
and an asymmetrical shape will result in re-radiation 
of absorbed sunlight that is uneven between the two 
hemispheres. This disparity results in a net acceleration 
(or deceleration) of the satellite's orbit and can there- 
fore cause orbital evolution. This effect has not been 
observationally verified, and in this section we evalu- 
ate BYORP 's relevance in explaining the observed spin- 
orbital characteristics of NEA systems by introducing 
current theoretical models predicting BYORP 's effects 
and timescales. 

For all observed asynchronous satellites, we rule out 
BYORP as a major player in their recent evolution be- 
cause this effect depends on spin-orbit synchronization; 
without synchronous spin, the effect cancels out. This is 
applicable to nearly half of the satellites in our sample 
that are asynchronous (Tabled]): 2004 DC, 2003 YTl, 
1991 VH, 2001 SN263 #2, and 1994 CC #2. We 
now examine scenarios in which BYORP may have had 
an important role in the past evolutionary histories of 
these systems, assuming their satellites used to be syn- 
chronous. 

First, we consider the case where a synchronous satel- 
lite evolved via BYORP and broke spin-lock through 
a planetary encounter. Such a scenario would imply 
that its past BYORP-affected orbital properties would 
be erased in part by the scattering encounter; therefore, 
any observed properties cannot be wholly attributed to 
BYORP, but would also be attributed to the flyby. 

Second, we consider the case where BYORP expands 
a synchronous satellite's orbit and increases the eccen- 
tricity enough to break spin-lock and result in chaotic 
asynchronous rotation (jCuk fc Nesvorn-^ |2010[ ). This 
could perhaps explain the observed eccentric and asyn- 
chronous satellites, but we find this unlikely for two rea- 
sons. First, synchronous re-capture is thought to oc- 
cur rapidly (jCuk fc Nesvorn-^120101 on the order of ^10^ 
years) and would prevent a substantial population of 
asynchronous binaries from forming. Second, this pro- 
cess would not explain how all observed asynchronous 
satellites acquired spin periods much less than their re- 
spective orbital periods. 

In essence, although we cannot rule out that BYORP 
played a role in the evolution of asynchronous satellites, 
we find that tides (for spin synchronization) or plane- 
tary encounters are required to explain the data, whereas 
BYORP is not. Therefore, we find that BYORP alone 
cannot explain the properties of asynchronous satellites, 
which includes all satellites with large semi-major axes. 
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We apply the same reasoning to widely-separated sys- 
tems such as 1998 ST27 and the outer satellites of the 
triples. Because neither tides alone nor BYORP alone 
can readily explain the origin of large semi-major axes, 
we conclude that these properties may be primarily a 
result of post-fission dynamics rather than evolutionary 
processes. 

For all remaining synchronous satellites, BYORP may 
be responsible for their observed spin-orbital characteris- 
tics. However, BYORP is neither necessary nor sufficient 
to explain these properties: another evolutionary pro- 
cess (tides) is required to first synchronize the satellite's 
spin before BYORP can operate. BYORP 's relevance for 
NEA systems is also complicated by the short BYORP 
timescales. An important unresolved issue with our un- 
derstanding of BYORP has to do with confiicting models 
on the sign of de/dt with respect to da/dt. These issues 
are discussed in the following paragraphs. 

Current theories predict different behaviors of semi- 
major axis and eccentricity evolution due to BY- 
ORP (IGoldreich fc Saril [20091 : ICuk fc Nesvornyl [20T0I : 
McMahon fc Scheer cs '2010a":c': ' Steinberg fcTanr[20T¥ . 
McMahon fc Scheeres (,2010a ) describe how the overall 
orbital evolution over long timescales due to BYORP 
causes the semi-major axis and eccentricity to evolve 
in opposite directions, even with the inclusion of libra- 
tion effects. ICuk fc Nesvorn-s^ (|2010l ) also include the 
effects of a secondary's librations due to its elongated 
shape and predict that librations dominate over direct 
perturbations by BYORP. They describe how the semi- 
major axis and eccentricity evolve in the same direc- 
tion. They suggest that the most likely outcome of an 
initially expanding orbit (with e > in their model) 
is chaotic rotation of the secondary followed by syn- 
chronous spin re-establishment and inward migration, 
which would prevent ey olution to large semi-major axes. 
[Steinberg fc Saril (|2011f ) also find that the preferred end 
state for BYORP is shrinkage of the sem i -major axis. 
On the other hand, IMcMahon fc Scheered (|2010bf ) sug- 
gest that librational moti on for binaries like 1999 KW4 
(~5° amplitude librations: lOstro et al.ll2006l ) will remain 
small for expanding orbits (with e < in their model), 
preventing chaotic rotation of the secondary, and leading 
to large semi-major axes. Future measurements of spin 
states of satellites can help constrain which model best 
captures the correct behavior. 

The timescales for BYORP evolution are thought to 
be fast. ICuk fc Burns! (|2005f ) find that for most asteroid 
shapes, the BYORP torque is significant and can modify 
the satellite's orbital semi-major axis and eccentricity on 
timescales of ^^10^ years. Similar expansion time scales 
have been found bv IMcMahon fc Scheered (|2010a[) . who 
assert that orbits can typically expand to their Hill 
sphere due to BYORP on the order of 10^ — 10® years. 
Currently, 1999 KW4 is the only well-characterized NEA 
binary with published shapes of both pri mary and sec- 
ondary, and lMcMahon fc ScheeresI (l2010af ) predict an or- 
bit expa nsion rate of 7 cm yr~^ (a prediction corrobo- 
rated bv [Steinberg fc Saril (|2011[ )) and an orbit-doubling 
time of ^22,000 years. Under the influence of BYORP, 
primary oblatcness, and solar perturbations, they pre- 
dict that 1999 KW4 will reach its Hill radius in -54,000 
years. Incorporation of librational effects into the model 



by IMcMahon fc ScheeresI (|2010ar) causes a longer BY- 
ORP evolution, but the mutual orbit will still expand 
to it s Hill radius for sma ll librational motion to at least 
12° ([McMahon fc Scheer cs 2010b). 

There are potential puzzles with the hypothesized 
short BYORP timescales and assertions that BYORP- 
induced orbit expansion can reach the Hill radius. First, 
binary formation would need to be rapid enough to 
produce the observed fraction of NEAs with satellites. 
Since NEAs have dy namical lif etimes on the order of 
a few million years ([Bottke et a l. 2002), much larger 
than the BYORP disruption timescale of '-^O.l million 
years, the binary production rate would need to match 
the rate of binaries disrupted by BYORP for a steady- 
state binary NEA population. Part of this issue (at 
least for nominally half of all systems in which BY- 
ORP causes inward migration) may be mitigated if a 
stable equilibrium exists bet ween BYORP and tides (e.g 
iJacobson fc Scheere£|[20lTbl) . If BYORP acts to contract 
the orbit and tides cause expansion of the orbit, there 
is a critical semi-major axis at which their effects may 
balance and result in a stable equilibrium. Second, ac- 
cording to McMahon & Scheercs (2010a), BYORP can 
expand the orbits (nominally half) of satellites to their 
Hill radii. This means we should observe some systems 
that are at significant fractions of their Hill radii; for 
instance, the Hill radii of these systems are typically 
a few hundred primary radii. However, we do not ob- 
serve any systems wider than 20 primary radii. More- 
over, such BYORP-induced expansion to the Hill sphere 
will cause the expanding binary to become more sus- 
ceptible to planetary encounters by shortening the en- 
counter timescale and strengthening the perturbative ef- 
fects of enc ounters. Using equations and timescales in 
[Fang fc Mar got (20n|), we calculate that if 1999 KW4 
reaches a separation halfway to its Hill radius. Earth en- 
counters at a typical encounter velocity of 12 km are 
expected to disrupt the binary at encounter distances of 
—28 Earth radii, which can occur every —260,000 years. 
With a separation of 90% of its Hill radius, disruptive 
encounters can occur every —100,000 years. As a result, 
we suggest that planetary encounters may also create as- 
teroid pairs in the near-Earth population. 

We also mention that the predictions for 
1999 KW4's expansi o n to i ts Hill radius given by 
IMcMahon fc ScheeresI ([2010af) are not ruled out by 
our calculations of upper limit critical semi-major axes 
ac,peim at which a synchronous satellite will have its spin 
lock broken. We calculated these critical semi-major 
axes (Table [3]) using Equation [71 which compares the 
torque due to a satellite's permanent deformation 
(which maintains synchronization in the absence of 
significant eccentricities) and the torque due to YORP 
(which can spin an asteroid away from synchronization) . 
Table O shows that for all NEA satellites, the critical 
semi-major axis Ocporm is much larger than the Hill 
radius, meaning that synchronization can be maintained 
in the absenc e of other perturbations. T hus, if the 
predictions by IMcMahon fc Scheered (I2010a[) are correct 
(i.e. eccentricities remain small when the orbit expands), 
the formation of asteroid pairs by BYORP-induced 
orbit expansion (and possibly planetary encounters) is 
possible. 

To summarize, we find that BYORP is not currently 
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acting for nearly half of the satellites in our sample 
of NEA binaries and triples because they are asyn- 
chronously rotating. We also find difficulties associated 
with the hypothesis that BYORP operated in the past 
evolutionary histories of these currently asynchronous 
satellites, although we cannot rule out that their orbital 
properties may have been shaped in part by BYORP. 
For the remaining synchronous satellites, BYORP may 
be acting but is not required to explain observed orbital 
properties. 

4. EVOLUTION BY CLOSE PLANETARY ENCOUNTERS 

In this section, we examine if close planetary encoun- 
ters can explain the observed spin-orbital properties of 
well-characterized NEA systems. 

The majority of satellites in our sample of NEA sys- 
tems are synchronously rotating and on near-circular or- 
bits. This includes 2000 DP107, 1999 KW4, 2002 CE26, 
Didymos, and the inner satellites of 2001 SN263 and 
1994 CC. We find that close planetary encounters are un- 
likely to explain these properties. We expect that plan- 
etary encounters are strong enough to cau se changes to 
an asteroid's spin state fe.g. lScheeres et al. 2000 20041 
and orbital eccentricity (e.g. iFang fc Margotll2011l ). such 
that we would not expect a majority of synchronous and 
circular binaries if planetary flybys were a dominant pro- 
cess. Similarly, planetary flybys are also not likely to ex- 
plain the spin-orbital state of 2001 SN263's outer satel- 
lite, which although asynchronous, has a circular orbit. 

For asynchronous binaries 2003 YTl, 2004 DC, and 
1991 VH, their eccentricities are 0.18, 0.3, and 0.06, re- 
spectively. The outer satellite of 1994 CC is also asyn- 
chronous with an orbital eccentricity of 0.19. 2003 YTl, 
2004 DC, and 1991 VH have semi-major axes that in- 
dicate their expected eccentricities following post-fission 
dynamics were higher or comparable to their current 
eccentricities (Figure [5]). If they formed at their cur- 
rent semi-major axes, either tides did not have sufficient 
time to damp the eccentricities to circular orbits (Ta- 
ble [2]) , or a planetary encounter may have erased some 
of the tidal damping effects by increasing the eccentrici- 
ties. If they formed closcr-in (with correspondingly lower 
eccentricities following post-fission dynamics), a plane- 
tary encounter may have increased both their semi-major 
axis and tidally-damped eccentricity to observed values. 
1994 CC's outer satellite has tidal damping timescales 
much longer than its total possible lifetime, and so its 
current eccentricity is too low to be explained by post- 
fission dynamics and tides. In this case, it is possible 
that planetary encounters may have lowered its eccen- 
tricity from a high predicted value following post-fission 
dynamics to its observed value. 

To investigate these scenarios, we employ results from 
a companion paper on the effect of plaiietary encounters 
with binary asteroids ()Fang fc Margotir2011|) . We find 
that for a planetary flyby to increase each of the asyn- 
chronous binaries' eccentricities from a tidally-damped 
value of to its observed value would take ^4.94 Myr 
for 2003 YTl, -6.54 Myr for 2004 DC, and -0.56 Myr 
for 1991 VH. These encounter timescales represent the 
average time for an eccentricity increase when the binary 
is in near-Earth space, assuming the NEA follows a tra- 
jectory from main belt source regions to its current orbit 
in near-Earth space. These timescales can occur within 



the near-Earth dynamical lifetime on the order of a few 
million years. Inclusion of additional planets and repeat 
passes will make it more likely th at planetary encounters 
can affect orbital properties (see IFang k, MargotI (|2011[ ) 
for details). For the case of 1994 CC with a starting ec- 
centricity of 0.8, the largest decreases in eccentricity hap- 
pen for shorter encounter distances. However, at shorter 
encounter distances, unstable encounters with collisions 
and ejections dominate. For instance, at an encounter 
distance of 2 Earth radii, the average post-encounter ec- 
centricity is —0.58 (not low enough to match the observed 
eccentricity of 0.19) and the percentage of stable encoun- 
ters is less than 1%. Smaller encounter distances are even 
more problematic. Scenarios with repeat passes, where 
the eccentricity has a net decrease to its observed value, 
are also very unlikely. Therefore, our results suggest that 
planetary encounters cannot decrease the orbital eccen- 
tricity of 1994 CC's outer satellite from high values to 
moderate values. 

Next, we investigate if the widest orbits observed in 
NEA systems originate from post-fission dynamics (Sec- 
tion 11.31 and Figure or possibly another mechanism, 
namely planetary encounters. We examine a hypotheti- 
cal binary modeled after 1998 ST27 with a large separa- 
tion of —16 primary radii and an eccentricity of 0.3. Our 
simulations used encounter distances of 2—12 and 
encounter velocities of 8—24 km s~^ for an Earth-mass 
perturber. We find that encounters cannot create wide 
systems, except very rarely. Here are a few illustrative 
cases we examine: 

1. We perform simulations starting with a circular bi- 
nary with a separation of 4 primary radii, which is typical 
for NEA binaries. For a typical encounter velocity of 12 
km s~^, encounters at a distance of 8 i?© showed that 
—7% of stable encounters at least doubled in semi-major 
axis and none of the stable binaries quadrupled in semi- 
major axis. Encounters resulting in stable binaries that 
at least doubled in semi-major axis had post-fiyby ec- 
centricities of > 0.45, which is larger than any observed 
eccentricity but may damp within a dynamical lifetime. 

2. Wc perform a "repeat encounter" scenario starting 
with an eccentric binary (e = 0.5, 0.7) with a separation 
of 8 primary radii. These initial conditions assume that 
an a prior encounter already doubled the semi-major axis 
from 4 to 8 primary radii and increased the eccentricity to 
high values. We find that it is possible to both decrease 
the eccentricity and further increase the semi-major axis 
to observed values. However, this only occurred in —5% 
of stable encounters with a fine-tuned initial semi-major 
axis and eccentricity in order to match the observed val- 
ues. 

3. We investigate if this hypothetical wide binary used 
to be a triple system by performing simulations of a pri- 
mary asteroid and 2 satellites initially located at 4 and 16 
primary radii with circular orbits. The chance of eject- 
ing the inner satellite with the outer satellite intact is 
possible, but only occurs in <3% of all unstable encoun- 
ters with encounter distances less than — 8i?0. When it 
does happen, the final eccentricity of the outer satellite 
is typically higher than any observed eccentricity. 

Our work shows that planetary encounters are not nec- 
essary to explain the observed spin-orbital properties of 
binary NEAs, but in reality, close planetary encounters 
will occur and they will change the orbital characteris- 
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tics of NEA multiples (|Fang fc Margotll201lD . To recon- 
cile these two results, there are a few possibilities: (a) 
the observed synchronous and circular systems are more 
recent migrants to ncar-Earth space that have not yet un- 
dergone deep planetary encounters, and oppositely, the 
asynchronous and eccentric systems have been interact- 
ing with planets for longer periods of time, or (b) the 
observed systems have encountered terrestrial planets, 
but tidal evolution occurs on a faster timescale, and has 
managed to synchronize and circularize some systems. 
Possibility (b) may or may not be supported by calcula- 
tions of tidal damping timescales (Table [2]) depending on 
the tidal Love number model and NEA-spccific encounter 
timescales. 

In summary, we find that planetary encounters cannot 
explain observed synchronous satellites nor nearly circu- 
lar orbits. For asynchronous satellites with non-circular 
orbits (2004 DC, 2003 YTl, and 1991 VH), we find that 
planetary encounters are a plausible explanation that can 
reproduce their observed eccentricities, although not nec- 
essary to explain the data (observed properties can be 
explained by tides alone). However, planetary encoun- 
ters cannot decrease 1994 CC's eccentricity from a high 
post-fission dynamics value to the observed value. We 
also find that the observed wide orbits were unlikely to 
be reached through planetary encounters. 

5. EVOLUTION BY KOZAI RESONANCE 

Now we examin e orbital evolution by the Kozai reso- 
nance (|Kozaill962( ) by reviewing its effects and timescales 
as well as its appli cability to o bserved NEA systems. 

Kozai resonance (|Kozailll962D is a secular effect causing 
angular momentum exchange between an inner body's 
orbital eccentricity and relative inclination with a mas- 
sive, outer perturbcr. In this resonance, eccentric- 
ity e and inclination / are coupled, and the quantity 
Vl — cos/ is conserved in the idealistic case of an 
outer perturber with a circular orbit. For an initially 
circular binary, large Kozai oscillations can occur if the 
relative inclination between the inner and outer orbits is 
at least ^39.2 degrees, and this critical inclination de- 
pends on the ratio of their semi-major axes. Since this 
mechanism can result in high eccentricities, we consider 
Kozai resonance betwe en an NEA binary and the Sun 
as the outer perturber (jPerets fc Naozll2009l) . Trip le sys- 
tems are not considered here, as iFang et al.l ()201lD have 
shown that the presence of a second satellite damps any 
Kozai oscillations in the system and thus will not be a 
relevant mechanism in modifying the satellites' eccentric- 
ities. 

For binaries under the influence of the Kozai mecha- 
nism, we can estimate the Kozai period Pxczai or typical 
osc illation timescale be tween limiting values of e and / 
as (jKiseleva et al.lll998[) 
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where Pq is the heliocentric orbital period, i-l,inary is the 
binary's mutual orbital period, and Bq is the heliocentric 
eccentricity. The binary's primary mass is Mp, the sec- 
ondary mass is M^, and the outer perturber's mass is the 
Sun's mass {Mq) in the situation we examine here. In 
the absence of other effects, we calculate the Kozai peri- 



ods for all well-characterized NEA binaries and present 
them in Table |4l Binaries with small heliocentric dis- 
tances and/or large heliocentric eccentricities (such as 
1999 KW4) have short oscillation periods. The very short 
timescales in Table|4]indicate that evolution by the Kozai 
mechanism is faster than any other evolutionary process 
examined in this study. 

Determination of whether binaries are in the Kozai 
regime or not can be assessed for at least 2 near-circular 
binaries in our sample, 2000 DP107 and 1999 KW4, for 
which we have detailed information about the orientation 
of the mutual orbital plane. For these binaries, we calcu- 
late the current inclination between the binary's mutual 
orbit and the Sun's apparent orbit. We find that neither 
1999 KW4 nor 2000 DP107 have inclinations that meet 
the critical Kozai angle. Other binaries in our sample 
do not currently have reliable orbital orientations. The 
difficulty of measuring precise orbital plane orientations 
without sufficient radar coverage makes the assessment 
of Kozai infiuences difficult to verify at the moment, and 
additional observations of NEA binaries are encouraged. 

Fulfillment of the required Kozai inclination will be 
affected by processes that can change the relative incli- 
nation between the binary's mutual orbit and its helio- 
centric orbit. For a satellite in an equatorial orbit with 
a sem i-major axis of several primary radii. ICuk fc Burnsl 
(j2005| ) describe a possible equilibrium state where YORP 
and BYORP torques balance. They describe how this 
scenario can occur when the inclination between the 
primary's equatorial plane and the heliocentric orbit is 
~50— 60 degrees. If the binary's mutual orbit is in the 
same plane as the primary's equator (as would be ex- 
pected from the generally accepted rotational fission for- 
mation model), then the stable inclination of ^50—60 
degrees would be the angle between the bi nary orbit and 
heliocentric orbit. [Steinberg fc Saril ()2011[ ) describe a dif- 
ferent stable inclination that will be reached under the 
effects of BYORP, which will tend to orient a binary's 
orbit into an inclination of either or 90 degrees relative 
to its heliocentric orbit. However, we note that in the ab- 
sence of Kozai-damping processes, a "stable" inclination 
of ^50—60 or 90 degrees would cause a binary to undergo 
Kozai cycles. Although the value of the stable inclination 
is unc lear due to d ifferent predictions by ICuk fc Burnsl 
(|2005l) and lSteinberg k Sari (,20llT ). we suggest that ra- 
diative torques may cause a binary that is initially not 
affected by Kozai effects to enter the Kozai regime when 
its relative inclination with the Sun's apparent orbit is 
sufficiently high. 

Next, we consider the effect of primary oblateness (de- 
scribed by a J2 term) in modulating the Kozai effect 
for NEA binaries. The critical semi-major axis ac,j2 for 
the transition between the influence regions of primary 
oblatenes s and solar dynamics such as the Kozai reso- 
nance is (jNicholson et al.|[2008l ) 
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where J2 approximates the non-spherical shape of the 
primary by its level of oblateness, Mp is the mass of the 
primary, Mq is the mass of the Sun, Rp is the primary's 
radius, and Qq is the heliocentric semi-major axis. We 
calculate ac,j2 for 3-11 NEA binaries in our sample for J2 
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Table 4 

Kozai Oscillation Periods 



Binary 


a (km) 


e 


-PKozai (yr) 


2000 DP107 


2.62 


0.01 


89.11 


1999 KW4 


2.55 


0.008 


10.78 


2002 CE26 


4.87 


0.025 


756.60 


*2004 DC 


0.75 


0.30 


269.78 


*2003 YTl 


3.93 


0.18 


58.40 


Didymos 


1.18 


0.04 


545.99 


*1991 VH 


3.26 


0.06 


81.34 



For each binary, adopted values for the ob- 
served semi-major axis a and eccentricity e 
as well as the Kozai oscillation period PKozai 
are listed. 

* Known asynchronous satellites 

values ranging from 0.001 to 0.1, and we find that the 
range of possible Ccjj values is larger than the observed 
semi-major axis for all binaries. This implies that these 
binaries are dominated by primary oblateness, which will 
cause orbital precession that can completely suppress 
Kozai cycles. 

Lastly, we note that this mechanism can alter an NEA 
binary's spin state and eccentricity, but cannot change 
the semi-major axis of the binary's mutual orbit in the 
averaged problem. Therefore, it cannot explain the ob- 
served wide orbits, which again suggests that these large 
separations arc inherited from post-fission dynamics. 

To summarize, Kozai resonance is not applicable for 
systems where its effects may be damped, and this in- 
cludes all satellites in triple systems and binaries where 
primary oblateness may dominate. This includes all sys- 
tems in our sample of well-characterized binaries and 
triples. The averaged Kozai effect also cannot modify 
the semi-major axis. 

6. COMBINED ORIGIN AND EVOLUTION OF 
SPIN-ORBITAL PROPERTIES 

Here we summarize the relevance of each evolutionary 
process examined in this study (tides, BYORP, planetary 
encounters, and Kozai effects) in the context of explain- 
ing the origin of observed spin-orbital properties of NEA 
systems in our sample (Table [T]). Then we discuss evo- 
lutionary pathways between observed spin-orbital states 
and disrupted binaries (asteroid pairs and contact bina- 
ries). 

Table [5] shows how each considered evolutionary pro- 
cess matches up against each other in explaining the spin- 
orbital origin of binaries and triples in our sample. Tidal 
evolution is a dominant process that can explain the 
satellite's spin state and orbital eccentricity for nearly all 
systems assuming formation by rotational fission followed 
by post-fission dynamics. All other examined processes 
either depend on tidal evolution and/or are not required 
to explain the observed systems. For instance, BYORP 
requires a tidally-synchronized system to operate. Al- 
though BYORP is not required to explain the data, it 
may be acting in synchronous systems or may be pivotal 
in a stable equilibriu m state where the effects of tid es and 
BYORP cancel (i.e. iJacobson fc ScheeresI IMlB . Fly- 
byS; or planetary encounters, can potentially explain sev- 
eral asynchronous and non-circular systems, but not the 
majority of circular and synchronous systems. The Kozai 



effect is not applicable for primary oblateness-dominated 
NEA systems in our sample. We also find that ob- 
served wide orbits, such as those of binary 1998 ST27 
and the outer satellites of 2001 SN263 and 1994 CC, 
are most likely a dire ct byproduct of post-fission dynam- 
ics (jJacobson fc Schccrcs.201lT ) because none of the four 
evolutionary processes as currently modeled seem capa- 
ble of delivering satellites to such large separations. 

Now we discuss these evolutionary processes as an en- 
semble of pathways between spin-orbital states and dis- 
rupted states such as asteroid pairs and contact binaries. 
These pathways are shown in Figure [3] with the assump- 
tion that these systems are not affected by the Kozai 
effect. In the figure, the yellow box is the starting point 
for a rotationally fissioning asteroid. Green boxes repre- 
sent potential outcomes. The evolutionary pathways are 
drawn in solid black (for post- fission dynamics), solid 
gray (for planetary encounters), and dotted black (for 
BYORP), with red lines representing differences between 
the tw o different BYORP models. iCuk fc Nesvorn\^ 
(|2010[ ) predict that BYORP-induced change in eccentric- 
ity has the same direction as the semi-major axis change 
(i.e. when t he semi-major axis inc r eases, the eccentricity 
increases). iMcMahon fc Scheered ()2010al O) predict that 
the eccentricity moves in the direction opposite to the 
semi-major axis change (i.e. when the semi-major axis 
increases, the eccentricity decreases). 

As adopted in Figure [31 the term asynchronous can 
refer to two types of states where the satellite's spin is 
non-resonant with its orbital period: "quasiperiodic" ro- 
tation, which is stable, and "chaotic" rotation, which 
is unstable. When these spin state trajectories are de- 
picted on surface of section plots, "quasiperiodic" trajec- 
tories are smooth curves and "chaotic" trajectories are 
random points filling in an area on the plot. The sta- 
ble or unstable nature of an asynchronous rotation can 
be readily di scerned on such surface of section plots, as 
discussed in iMurrav fc DermofrO (Il999f) and shown for 
Saturn's satellite Hvperion in iWisdom et all ()1984f ). A 
satellite's principal moments of inertia {B — A)/C deter- 
mine the strength of spin-orbit resonances. The higher 
the eccentricity, the lower the necessary (B — A)/C at 
which chaotic rotation can occur. An increase in eccen- 
tricity (due to perturbations such as BYORP) can lead 
to overlap between spin-orbit resonances. Simultaneous 
libration in more than one spin-orbit resonance is not 
possible, and this leads to chaotic behavior. A satellite 
can depart an unstable, "chaotic" state by reaching a 
stable "quasiperiodic" state and temporarily re-establish 
synchronization. In some models, this allows BYORP to 
regain control of the system and evolve it away from the 
chaotic states, as shown in simulations for 1999 KW4-like 
binaries bv lCuk fc Nesvornvl (|2010f l. 

In Figure [31 we discuss spin-orbital pathways that 
are common to both BYORP models (lines not shown 
in red). Post-fission dynamics (<1000 years) combined 
with tidal dissipation (and possibly BYORP) can pro- 
duce many observed examples of NEA systems, including 
both synchronous and asynchronous rotators, both cir- 
cular and eccentric orbits, contact binaries, and asteroid 
pairs. Evolution between spin-orbital outcomes (high- 
lighted in green in the figure) is possible through close 
planetary encounters; this can cause both asynchronous 
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Table 5 

Origin of Spin-Orbital Properties 



System 


a/Rp 


e 


Tides 


BYORP 


Flyby 


Kozai 


2000 DP107 


6.6 


0.01 


Yes 


Yes 


No 


No 


1999 KW4 


3.9 


0.008 


Yes 


Yes 


No 


No 


2002 CE26 


2.8 


0.025 


Yes 


Yes 


No 


No 


*2004 DC 


4.4 


0.30 


Yes 


No 


Yes 


No 


*2003 YTl 


7.1 


0.18 


Yes 


No 


Yes 


No 


Didymos 


3.0 


0.04 


Yes 


Yes 


No 


No 


*1991 VH 


5.4 


0.06 


Yes 


No 


Yes 


No 


2001 SN263 #1 


2.9 


0.016 


Yes 


Yes 


No 


No 


*2001 SN263 #2 


13 


0.015 


Yes 


No 


No 


No 


1994 CC #1 


5.6 


0.002 


Yes 


Yes 


No 


No 


*1994 CC #2 


20 


0.19 


No 


No 


No 


No 


For each system's 


satellite, 


we list the adopted values for the observed 



semi-major axis a/ Ftp and eccentricity e, as well as which evolutionary 
process(es) can explain its observed spin-orbital state. 'Yes' (in green) 
means that the considered process is plausible and we find no evidence 
to rule it out; the opposite response is 'No,' (in red) meaning the 
considered process is highly unlikely. 
* Known asynchronous satellites 

and synchronous rotators in circular orbits to evolve into 
asynchronous satellites with eccentric orbits, or any state 
into an asteroid pair or contact binary (shown by arrows 
with no starting point in the figure). 

We now consider inward evolution from "synchronous 
and circular" and "asynchronous (chaotic) and eccentric" 
states. A synchronous satellite will be affected by BY- 
ORP and depending on its shape and rotation, can mi- 
grate inwards and become a contact binary. In the model 
of iCuk fc Nesvorn?! (|20100 . a chaotically-spinning (and 
therefore asynchronous) satellite in an eccentric orbit will 
eventually temporarily regain synchronization with an 
orientation switched by 180 degrees, in which case BY- 
ORP would regain control of the system's evolution and 
restore it to a stable state. Inward migration can fol- 
low and lead to the formation of a contact binary or, 
after sufficient inward migration, the effects of BYORP 
and tides can cancel, leading to a stable equilibrium (e.g. 
iJacobson fc Scheeresll2011br l and a tidally-damped eccen- 
tricity. This would result in a synchronous satellite in a 
circular orbit. 

We now consider outward evolution from the "syn- 
chronous and circular" and "asynchronous (chaotic) and 
eccentric" states. According to lCuk fc Nesvornvl (|2010f ). 
outward migration leads to an increase in eccentricity, 
so a synchronous and circular state can become chaoti- 
cally asynchronous and eccentric. Once chaotically asyn- 
chronous and eccentric, it would remain so unless the 
satellite can re-orient itself such that if synchronization is 
temporarily regained, BYORP can migrate the satell ite 
inwards. According to iMcMahon &: ScheeresI ()2010al O). 
outward migration leads to a decrease in eccentricity, 
which means that synchronization will be maintained, 
aided in part by the restoring torque due to the satel- 
lite's permanent deformation. In this model, outward 
migration would allow the "synchronous and circular" 
state to evolve into an asteroid pair. 

Examples of NEA satellites in each spin-orbital state 
are also given in Figure [31 There are several examples 
of each end state, except for asynchronous satellites with 



circular orbits and asteroid pairs. The only known ex- 
ample of an asynchronous satellite with a circular orbit 
is the outer satellite in triple 2001 SN263. Asynchronous 
and circular configurations may be rare because asyn- 
chronous satellites are more likely to be well-separated 
from their primary, and larger separations are more 
susceptible to stronger disruptiv e planetary e ncounters 
with short encounter timescales (jPang fc Marg ot 2011). 
Strong and frequent planetary encounters can easily 
evolve asynchronous and circular configurations to asyn- 
chronous and eccentric states, making it more rare to 
observe an asynchronous and circular system. 

As for asteroid pairs, none have been definitively 
reported in the near-Earth population a lthough they 
have been observed in the main belt (i.e. iPravec et al.l 
|2010[) . We briefly discuss implications for the forma- 
tion of NEA asteroid pairs. Figure [3] shows how as- 
teroid pairs can form directly through post-fission dy- 
namics or planetary encounter s. The BYORP model by 
IMcMahon fc Scheered (|2010al O) provides an additional 
route to their creation via BYORP-dominated expan- 
sion of a satellite's orbit and the BYORP model by 
iCuk fc Nesvorn-sl ()2010D does not predict this. Determi- 
nation of likely asteroid pair formation mechanisms can 
el ucidate which BYORP model (ICuk fc Nesvornv (2010|) 



or IMcMahon fc Scheerej ()2010al lcf)) predicts the correct 



behavior. If BYORP-induced expansion of a binary's 
mutual orbit is responsible for creating NEA asteroid 
pairs, then we should expect to observe a few binaries 
with very large separations unless they are efficiently dis- 
rupted by planetary flybys (see Section Alternatively, 

ICuk fc Nesvorn-sl (|2010| ) have hypothesized that asteroid 
pairs are the result of scattering among satellites in triple 
systems. Such scattering can occur if both satellites are 
or used to be synchronous and BYORP led them into 
convergent orbits, which would result in unstable orbits 
and scattering. Components would then collide or be 
ejected. 

7. CONCLUSION 
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Synchronous and circular 



Ex. 2000 DP107, 1999 KW4, 
1994 CC #1, 2001 SN263 #1 



McMahon & Scheeres (20101 
Outward migration by tides & BYORP 
Eccentricity damps by tides & BYORP 



Inward migration 
by BYORP 



Post-fission dynamics 

and tides 
(and maybe BYORP) 



Post-fission 
dynamics 




/ 



Planetary 
encounter 



Contact binary 
Ex. Castalia and Mithra 



Cuk & Nesvornv (20101 
Outward migration by tides & BYORP *^ 

Eccentricity grows by BYORP Cuk\ Nesvornv (20101 

Cliaotic librations breal< spin locl< chaotic sateiiite iibrations 

^ Sateiiite reorients & synchronizes 

s Inward migration by BYORP 

Post-fission ^ ^ 

dynamics «^ 



Ex. 2004 DC (?), 2003 YTl (?), 
1991 VH (?), 1994 CC #2 (?) 



Asynchronous ("chaotic") 
and eccentric 



Ex. 2004 DC (?), 2003 YTl (?), 
1991 VH (?), 1994 CC #2 (?) 



Figure 3. Possible spin-orbital pathways are shown for an evolving satellite in a newly post-fissioned NEA system 
with the assumption that Kozai oscillations are not relevant. Pathways that start and end in the same spin-orbital 
configuration are not shown in this figure, even though they may occur. Two different models are given for BYORP 
and are highlighted in red. iCuk fc Nesvorn-<1 (|2010l ) predict that BYORP-induced change in eccentricity has the same 
direction as the change in semi-major axis, and they model how a c haotica lly asynchronous satellite can reorient and 
regain synchronization. On the other hand, iMcMahon fc Scheered (l2010aB ) predict that BYORP-induced change in 
eccentricity has a direction opposite of the cha nge in semi-major axi s . For synchronous and circular systems, it is 
possible that they are in the stable equilibrium of lJacobson fc Scheeres! ()2011b[ ). See Section |6] for additional discussion 
regarding this figure. 



Radar observations of well-characterized NEA binaries 
and triples have uncovered a diverse set of spin-orbital 
properties: synchronous and asynchronously rotating 
satellites, circular and eccentric orbits, and a few large 
separations between the primary and the satellite. The 
formation of these satellites by rota tional fission followed 
by a post-fission dynamics model (jJacobson fc Scheeres! 
!2011h can produce asynchronous satellites with a vari- 
ety of primary separations and high orbital eccentricities. 
We investigated how a newly formed system can evolve 
to one of the observed systems by evaluating these evolu- 
tionary processes: tides, BYORP, planetary encounters, 
and Kozai effects. 

We found that post-fission dynamics and tides can ex- 
plain the observed semi-major axes, eccentricities, and 
satellite spin states of nearly all binaries and triples in 
our sample. Other evolutionary mechanisms do not ap- 
pear to be dominant processes: BYORP is not applica- 
ble to asynchronous systems and is not required to ex- 
plain the observed data, even though it may be acting; 
planetary encounters arc likely not responsible for cre- 
ating synchronous and circular configurations; and the 
Kozai effect will typically be suppressed by the primary's 
oblateness. 

We also illustrated the evolutionary pathways for satel- 
lites in binaries and triples after they have formed. 
Evolutionary processes such as tides, planetary encoun- 
ters, and BYORP can evolve a system between syn- 



chronous and circular, asynchronous and circular, and 
asynchronous and eccentric configurations. 
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